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The stress field at the tip of a crack of a thin plate of elastic material that is broken due to a 
mode III shear tearing has a universal form with a non-universal amplitude, known as the stress 
intensity factor, which depends on the crack length and the boundary conditions. We present in 
this paper exact analytic results for this stress intensity factor, thus enriching the small number of 
exact results that can be obtained within Linear Elastic Fracture Mechanics (LEFM). 



I. INTRODUCTION 

Linear Elastic Fracture Mechanics (LEFM) [l[ deals 
with the dynamics of cracks in elastic materials subject 
to the assumption that linear elasticity if applicable ev- 
erywhere in the material except for a very small region 
around the crack tip known as the process zone. While 
recent experiments and theories indicate that there exist 
examples where nonlinear elasticity may be important for 
a full characterization of the stress field before a dynami- 
cal crack, some exact results of LEFM continue to have a 
very important role in describing and understanding the 
dynamics of cracks in a variety of fracture contexts. In- 
deed, one of the most celebrated results of LEFM is the 
analytic description of the universal stress distribution in 
front of a crack. This analytic result is particularly sim- 
ple when the crack is developing under one of the three 
fundamental shearing modes that a material can suffer. 
The simplest form results from mode III cracking which 
is obtained under an out-of-plane shear tearing. In a 
cylindrical coordinate system r, 9, z, where r = coin- 
cide with the crack tip and 9 = ±7r with the crack faces, 
the stress in front of a crack takes on the form 
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where Km is the so-called "Stress Intensity Factor". 
The results (p} underlines the universality of the angu- 
lar dependence and of the square-root singularity of the 
stress at the crack tip 0, H[ . All the non- universal as- 
pects of the stress distribution are collected in the Stress 
Intensity Factor which depends on everything, including 
the crack length, the boundary conditions and the his- 
tory of the loads that drive the crack evolution. The 
value of the stress intensity factor determines whether 
the crack propagates or not; Irwin [3] established a rela- 
tionship between the energy release rate and the stress 
intensity factor for a crack growth. A crack propagates 
as long as the energy release rate into the crack tip bal- 
ance the dissipation involved in the crack propagation, 
i.e. the Griffith criterion [l|, Q. 

Solving for mode III is relatively easy since the funda- 
mental equation to be dealt with is the Laplace equation 
Q rather than the bi-Lapalace equation which is the rel- 
evant one for modes I and II. On the other hand in 3- 
dimensional materials mode III is unstable and attempts 



to break in this mode quickly turn to cracks developing 
under effective modes I or II. This is not the case for thin 
plates where mode III can be sustained in experiments as 
had been demonstrated recently in Ref. 0- The theory 
of thin plates is however more cumbersome, using again 
the bi-Laplace operator even for mode III cracks. Nev- 
ertheless the problem can be treated analytically [8l4lfl| 
with the result that the leading singularity at the crack 
tip is higher, going like r~ 3 / 2 . In this case, the shear 
force close to the crack tip can be written as 
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we assign tilde on Km to distinguish it from the known 
stress intensity factor for this mode, h is the thickness of 
the plate. 

We may mention that the Kirchhoff plate theory pre- 
dicts the stresses very accurately for small deflection and 
at some distance from the crack tip irjh > 1) (cf. 
and references therein). However, for large deflection or 
near the crack tip (rlh < 0.1), the singularity of the 
stress field weakens HH i an d therefore it may diverge 
from the actual results. 

The aim of the present paper is to present further an- 
alytic progress that culminates in the calculation of the 
stress intensity factor for cracks developing in thin sheets. 
The structure of the paper is as follows: in Sect. [TTJ we 
set up the problem for mathematical analysis. In Sect. 
Mil we study the properties of the analytic functions that 
need to be solved for and the boundary conditions. In 
Sect. IIVI we introduce the conformal map from the strip 
to the upper half plane which is used to advantage in 
simplifying the solution of the problem. Sect. [V] wraps 
up the problem and presents the analytic result for the 
Stress Intensity Factor. 



II. SETTING UP THE PROBLEM 

Consider a semi-infinite strip of a thin sheet with a cut 
in the in the middle of the left- most end of the strip, along 
the symmetry axis. The left boundary is subjected to 
anti-symmetric shear load fb, i.e. mode III configuration, 
cf. Fig. [TJ For small deflection, we can neglect the in- 
plane deformation of the middle surface. According to 
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FIG. 1: The mode III configuration. A semi-infinite strip of 
width 2a, and a cut (C-D) of length b in the middle of the 
strip. The upper left side of the strip (B-C) is subjected to an 
anti-plane positive load, and the lower left side to an identical 
but negative load. 

the Kirchhoff plate theory [TJj], the deflection w in the 
z-axis for pure bending becomes 



A 2 w = 
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where the operator 'delta' stands for the Laplace operator 
V • V; the bending moments and the twisting moment 
are 
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and the anti-plane shear stress 
Qx = / cr xz dz = 



-D-(A W ), 



Qy 



cr yz dz = -D—(Aw). 
dy y 
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where the sheet is on the xy plane, and the bending oc- 
curs along the z direction. Uij is a component of the 
stress tensor, D = Eh 3 /12(1 — v 2 ) is the so-called the 
flexural rigidity of a plate, h is the thickness of the plate, 
and E, v are the Young's modulus and Poisson's ratio, 
respectively. 

Introduce now the complex notation z = x + iy, z — 
x — iy, such that 



7!= = \{dx + id y ) = dz. 
dz 2 
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Since w is a solution of the bi-Laplace equation, w can 
be expressed as [TJ[ 



where 5R denotes "the real part", and (f(z), x( z ) are both 
analytic functions, accept maybe on the boundary. 
The complex first derivative becomes 

d x w + idyW = 2dzW 



= zy'{z)+x'{z) + v{z). (9) 
The bending moments and the twisting moment are 
M x = -D [(1 - ^)SR {zip" + X ") + 2(1 + v)$l (if/)] (10) 

My = -d [-(i- m (V+x'O + 2(1 + m (v'Wi) 

M xv = D(l - i/)9f (zip" + x") ■ (12) 
Summing over the two bending moments gives 



M x + My = -D(l + v)Aw 

= -4D(1 + v)U((p'). 



(13) 



and the complex representation of the shear force be- 
comes 



Qx + iQ y = -ADip"(z) (14) 
Next, the shear force exerted normal to the boundary 
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where n, s are the normal and the tangent direction ori- 
ented with respect to each other as the axes x, y. Multi- 
ply the last equation by ds and using the complex repre- 
sentation in Eq. (fH)) gives 



} n ds = Q x dy - Qydx 



2Di(ip"dz - tp"dz). 



(16) 



III. PROPERTIES OF <p AND BOUNDARY 
CONDITIONS 



In the case of mode III symmetry, the load and the 
deflection are anti-symmetric with respect to the x axis, 
(cf. Fig. [TJ) We consider only one crack lying on part 
of the x axis (C-D). The upper (B-C) and the lower left 
edges are displaced with a constant distance from each 
other, in a direction perpendicular to the strip plane. 
Then, w is anti-symmetric, but also Aw — 4^t[ip'(z)] is 
anti-symmetric. This implies that 



*%'(*)] = -»[¥>'(*)] 
the left and the right members give 



ip'(z) + ip'(z) = -ip'(z)-ip'(z) 
which may be arranged to 
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(18) 



w(z) =$t[zcp(z)+x(z)], 



(8) 



<p'(z) + y'(z) = -ip'(z) + ip'(z). 



(19) 
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Now, if f{z) is analytic on both the upper and the 
lower half plane, though not necessary on the boundary, 
the function (p(z) is also analytic on the same domain. 
Thus the left side of Eq. (|T9l is analytic in the whole 
domain (excluding the boundary). Since it equals to mi- 
nus its complex conjugate, it must be a pure imaginary 
function. However, according to the Cauchy- Riemann 
equations, an imaginary analytic function is a constant. 
Thus, 



<p'(z) = -<p'(z) + iC, 



(20) 



Where C is a real number. Taking the complex conjugate 
of this expression, followed by the complex conjugation 
of z, gives C = 0. Thus, the anti-symmetry with respect 
to the x axis prevails 



(21) 



Since <f'{z) is not necessary analytic on the x axis and 
on the boundary, it is advantageous to separate nota- 
tions, <p' + (z) and ip'_(z), for its parts on the upper and 
lower half-planes, respectively. 



A. Ahead of the crack 

In the case of a symmetric body with an anti- 
symmetric loading, it can be assumed that w(x, 0) = on 
the symmetry axis ahead of the crack (D-E). Moreover, 
the anti-symmetry of w dictates that Aw(x > 0,0) = 0. 
From Eqs. (O and ([2l]l 
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Aw + (a;,0) = 2Wff>'(z)} = <p'(x + i0) + tp'[x + i0) 

= ip (x + i0) + ip [x — iO) 
= <f'{x + iO) - <p'(x - iO) 
= ip' + (x) - ip'_(x) = (22) 

We used the subscripts + and — for the the limiting 
values of the functions reaching from the upper and the 
lower half of the plane, respectively. 

B. Crack face 

The crack face (C-D) is assumed to be traction free. 
Thus, the stress components acting normal to the crack's 
face become zero. This condition gives 

My = Mxy = Qy = 0, for C < x < D,y = (23) 

Integrating Eq. (|16p along the boundary from a point 
x on the crack face to end of the strip, i.e. point A — > oo, 
gives 

c-A c-A 

I Q n ds = 2Di I [ip" dz -7^dz\, (24) 
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Js = -2Di{ip'(x+iO)-ip'(x + iQ)-[ip'(oo)-Lp'(oo)]}. 

(25) 



Due to a finite slope at infinity, we can assume that 
lim^^oo if'{z = A) = 0. Also, on the crack face and on 
the upper free edge, Q y = 0. Thus, Eq. (f25j) becomes 



ip' {x + iQ) - ip' {x + iO) = — / Q n ds. (26) 



(27) 



and by pTj) . we find 

tp'+(x) + ip'_(x) =ij3, 

and 



2D 2D 



where ft, is the total applied force on the upper left 
boundary to create a constant deflection. 



C. The loaded edge 

This edge (B-C) is considered as a simply supported 
edge i.e. dw(—b, y)/dy = and M x (—b,y) = 0, for x 

h. y > 12]. These boundary conditions annul also 
M v (—b,y) — 0, and by the same consideration of Sec. 
(IIII Al) the complex function becomes 

<p' + (-b + iy)-ip'_(-b-iy) = (28) 

D. Free upper boundary 

As earlier, a free edge (B-A) feels a zero normal stress. 
Thus, 

My = M xy ()„ I) 
with the last condition giving 



ip (x + ia) — ip" (x + ia) = 0, for —b<x, y = a. (29) 



Applying the same integration as in Eq. (|25[) for the 
upper edge gives, 



(30) 



x+iO 



ip'(x + ia) — if'{x + ia) = 
and by Eq. (|2"Tj) . we obtain 

<f' + {x + ia) + <p'_(x — ia) = 0. (31) 
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FIG. 2: Mapping of the upper strip in Figfj] onto the upper 
half of £-plane. The inverse function is £ = cosh (^-^—^j, 
and the crack tip is at £o = cosh(7rfo/a) 



IV. CONFORMAL MAPPING 

The upper half of the strip is mapped onto the upper 
half of £, where £ = £ + it]. The mapping function is (cf. 
Fig. © 



z = Q(C) = -cosh _1 (0-* 
The analytic function can be written as 

and the boundary conditions become 



(32) 



(33) 



( F + (£)-F_(£)=0 
F + (£) + F_(£) = i/3 
F + (£)-F_(£)=0 

I f+(£) + f_(£) = o 



,r? = 



£d < £ 
1< £ <£d ,r? = 

-i< e < 1 ,r? = 

£ < -1 ,77 



In (34) 







These four equations constitute a Hilbert problem [13j . 
By introducing an auxiliary function G(£), we can write 
the boundary equations as a single equation. In our case, 
the function 



<?(£) = (C + i) 1/a (C-i) i/a (C-fr) 



1/2/ 



vl/2 



(35) 



with branch cut along 77 = and £ < — 1, 1 < £ < £d 
with a choice of a branch where G(£) — > £ 3 ^ 2 as £ — > 00. 
Then, 



G-(£) 



+1 for - 1< £ < 1, £ B < £ 



(36) 



G+(£) i -1 /or £<-l, 1<£<£d 
and equations (|34p can be written as a single equation 

where -ff (£) is Heaviside step function. 

Multiplication of the both sides with G+(£) gives 

G+(£)F + (£)-G_(£)F_(£)=iG + (£)/3[ J f7(£-l)- J f7(£-£ rj )] 

(38) 



Using the Plemelj's formulae [3j, we can define an an- 
alytic function in both upper and the lower half plane of 

c, 



■ D »G + (£)/3 



d£. 



(39) 



However, this solution of the Hilbert problem is not 
complete and we must add an analytic function in the 
whole plane, since it vanished in the left hand side of Eq. 
([57)1 . Thus, 



no 



2ttG(£) 



' D G + (£) 



d£ 



^o(C) 
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(40) 



where Po (£) is analytic in the whole plane and therefore, 
by Liouville's theorem a polynomial of finite degree. 

Since the moments are bounded as |£| — > 00, the degree 
of P(£) can not be greater than 1. In addition, upon the 
condition in Eq. (|2T|) the coefficients of the polynomial 
must be imaginary, i.e. Po(C) = *(PiC + Po)> where p\ 
and po are real numbers. 



V. SHEAR STRESS AND STRESS INTENSITY 
FACTOR 

From Eqs. (IT41 and (|3"3"|) the shear force in the £-plane 
becomes 



Q% 



F'(0 



-ADQip"(z)=ADQ-^, (41) 



and O'(C) = -— — is the map's derivative. 

Also, from the relation on Eq. (l2~Tj) , the shear force at 
the point i > £d along the £-axis becomes 



,(*) 



-2iD 



(42) 



here, we use the fact that fl'(t) = Cl'(t) in this section. 

Differentiate the solution of F(£) from the preceding 
section, Eq. (|40|) . gives 



*(PiC+Po)G'(C) , ipi 



(43) 



where 



/(C) 



G 2 (C) G(C) 

27TZ 7 X £-£ ? 



notice that G + (£) is imaginary. 

Now, the shear force away from the crack at the sym- 
metry axis is 
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G(t) 



(/3i(t)+pit+p o y(0i'(t)+pi) 

(44) 
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Applying the expression found in Eq. (1451) gives 



Qy(t) 



where 



4£>7T 



aVt - Ti 



G'(t) 
G(t) 



(45) 



G'(*) 



G(t) \t 2 -l 2(f-£c) 
I(t) is the Cauchy's principal value, 

and /'(£) is its derivative. 

The constants pi and po can be found from the bound- 
ary conditions. At the edge B — C, Q x (t) has the same 
expression as in Eq. (|4"5")) . Since we assume a finite ap- 
plied load at this edge, the singularity of the shear force 
at t — ±1 disappears if 



which gives 



pl(t)+p 1 t + p =0. 



Pi = 5 j8[I(-1)-/(1)], 
Po = -^[/(-l) +/(!)]. 



(46) 



(47) 



Finally, the stress intensity factor close to the crack tip 
in thin plates is obtained from 



K IH = Km(irx)*Q y (t), 
where x = Vt(t). For x 0, we can write 



Trsinh^ 

a 

thus, Eq. (l48l) becomes 



(48) 



(49) 



K 



in 



sinh 



7i b 



lim (t-f D )*Q y (t) 

t— HD 



(50) 



if/// = - 



2Z)7rVa 
sinh2(2i) 



(51) 



using (147)) 



if/// = 



- 1) 



zsinh5(f) Ji V^TV^I 



(52) 



notice that £|, — 1 = sinh 2 (^). Changing the variable 
of integration to A = ^ ^/ with the parameter fc 2 = 
f P 7i 1 one obtains 



sinh(^) f 1 



r/A 
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and finally, the stress intensity factor for mode III in thin 
sheets becomes 



Km 




sinh(^) 

7rb \ i 



cosh(H 



-K(k), 



(54) 



where 2a is the strip width, b is the crack length, /& is 
the total load and K(k) is the complete elliptic integral 
of the first kind. 

For a crack larger than the strip's width (b >> a), 
we can estimate the stress intensity factor in term of 
the deflection wq (^ << 1) instead of the load. In 
this case, K(k) ~ b/a and the stress intensity factor be- 
comes Km ~ fbb/y/a. An estimate of the load from 
the bending theory gives fb ~ D ^f a [H|, and we ob- 
tain Km ~ D ^ w ° ■ Hence, in particular, we see that 
the stress intensity factor decreases as the square of the 
crack length. 
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